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Abstract. Chaotic scattering with an internal degree of freedom and the possibility 
to generate directed transport of angular momentum is studied in a specific model, 
a magnetic dipole moving in a periodically modulated magnetic field confined to 
a compact region in space. We show that this system is an irregular scatterer in 
large parts of its parameter space. If in addition all spatio-temporal symmetries are 
broken, directed transport of mass as well as angular momentum occurs. The sensitive 
parameter dependence of the corresponding currents includes frequent sign reversals. 
Zeros of either quantity entail the exclusive occurrence of the other and thus give rise 
in particular to angular-momentum separation without mass transport as a classical 
analogue of spin-polarized currents. 



Chaotic angular-momentum pump 



2 



1. Introduction 

During two decades of research on chaotic scattering, internal degrees of freedom have 
not enjoyed much interest. Yet there are numerous reasons to give them a closer look: An 
internal freedom may provide the additional dimension required to render a scattering 
system with a single external degree of freedom chaotic [H El El E] . It may also act as a 
reservoir that absorbs energy from the external motion and stores it temporarily, thus 
giving the scattering process a transient inelastic character. Concerning applications, it 
is obvious that in chemical reactions — one of the paradigms of chaotic scattering [HIE]- 
the presence of internal freedoms is indispensable. In mesoscopic physics, in turn, the 
basic entities are electrons or fermionic quasi-particles: Their spin gives rise to a rich 
phenomenology beyond point-particle physics and presently receives much attention in 
the context of spintronics El E] . 

In electronic and other mesoscopic systems, transport is an issue of vital importance. 
It has been studied in depth in the framework of adiabatic pumps [IDlIHiEGlEin], electron 
counting [HI [15] and of the Kubo and Buttiker formulae [T6J [T71 [18]. In spintronics, 
the practical necessity to provide polarized currents has spurred research on directed 
spin transport El E]- In a different context, the study of ratchets [IH], inspired by 
the biological phenomenon of motor molecules [201 EU E2] , has revealed general physical 
conditions and numerous mechanisms for directed transport to occur in systems far 
from thermal equilibrium. Also here, internal freedoms are almost always involved 
and participate actively in transport processes [23], e.g., as functional (configurational) 
degrees of freedom in motor molecules [23, 25J and as internal heat baths. 

In this paper we present a case study of chaotic scattering with an internal degree 
of freedom, intrinsic angular momentum, and consider how it contributes to irregular 
scattering and how in turn directed transport of this quantity arises from scattering. The 
subject bears on various of the research lines mentioned above: We shall invoke much 
of what is known about chaotic scattering, in particular in periodically driven systems 
[26]. [27] . The coupling of a magnetic dipole to an inhomogeneous magnetic field [28] is 
crucial for our modelling. Also concerning chaotic dynamics in such a configuration we 
refer to earlier work [29]; in fact, it is the essential ingredient of one of the first models 
proposed for irregular scattering [30] . 

In directed transport, we take all the insights into account that concern the 
importance of binary spatio-temporal symmetries and their breaking [3T1 [32] [33] \M\ . 
We can here fall back on direct precursor work on chaotic pumps [35] . The definition of 
spin current is adopted from quantum contexts to define angular-momentum transport 
in a classical setting. It should be emphasized that we consider a driving that is 
both strong and fast [261 121] and thus not amenable to any perturbative and/or 
adiabatic approximation. In this sense, we are dealing with strongly nonlinear transport 
phenomena, far off the regime of linear response and in particular of peristaltic pumping 
[TOl [TT1 [T2l IT3] . As a consequence, directed transport is achieved already driving the 
system via a single parameter. The crucial role of chaos also distinguishes our work 
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from recent studies of spin ratchets in ballistic electron systems [31)1 I3"?l I3"8] . Preliminary 
results have been published on a national platform |39j . 

As our main achievement, we present conclusive evidence that in classical chaotic 
scattering systems, polarized currents can be generated, i.e., directed transport of 
angular momentum occurs independently of and even in the absence of mass (charge) 
transport. We explain this effect in terms of the interplay of asymmetry and the 
randomizing action of long unstable trajectories, thus corroborating the decisive role 
of chaotic scattering. It is also manifest in the sensitive parameter dependence of mass 
as well as angular-momentum currents, which allows for a fine tuning of both quantities. 
In particular, sign reversals are frequent and zeros of one current give rise to the exclusive 
dominance of the other. In this way, we achieve a "rectification of angular momentum" 
in the absence of mass transport. 

In El we construct our model and justify its peculiarities. We discuss details of the 
dipole-field coupling that have not been considered previously in the context of chaotic 
scattering. Section [3] is dedicated to irregular scattering as the basic dynamical category 
of the system. Various diagnostics are presented such as deflection functions, unstable 
periodic orbits, and time-delay statistics. Our central subject, directed transport, is 
approached in HI We define the relevant transport quantities in terms of reflection and 
transmission coefficients. Numerical evidence is provided for directed mass as well as 
for angular- momentum transport, and we point out the manifolds in parameter space 
where pure currents of either kind occur. Section [5] contains a summary and an outlook 
towards the quantization of our model with the perspective to construct a chaotic spin 
pump — evidently a major motivation for the present purely classical study. 

2. A model for chaotic scattering with intrinsic angular momentum 

2.1. Basic setup — coupling internal to external freedoms 

In constructing our model, we attempt to satisfy the following requirements, that (i) 
the system show chaotic scattering at least in parts of its phase space, as a dynamical 
mechanism for directed transport, (m) the potential couple the internal to the external 
degree of freedom and (Hi) vanish identically outside a compact scattering region so that 
numerical simulations can be restricted to a spatial box of finite extension, and (iv) be 
periodically time dependent in the form of delta kicks in order to facilitate reducing the 
dynamics to discrete time, that is, a map. 

Moreover, previous work on nonlinear transport mechanisms in ratchets and pumps 
has revealed a fundamental condition for directed currents to occur [31], [321 [331 EU 135] : 
an inhomogeneous phase space reflecting the absence of any spatio-temporal symmetry 
that would give rise to pairs of otherwise identical trajectories carrying the same current 
in opposite directions. It can be accomplished in particular if the dynamic is mixed 
(regular coexisting with chaotic motion) and the associated invariant manifolds exhibit 
some asymmetry. We take these additional criterion into account in setting up the 



Chaotic angular-momentum pump 



4 



model. 

Specifically, we choose the internal freedom as an intrinsic angular momentum (in 
the following we shall also refer to it as spin, wherever its classical nature is obvious). 
In connection with a charge it gives rise to a magnetic dipole moment that couples 
to an external magnetic field and thus, if the field is not homogeneous, to the spatial 
motion of the particle. However, we neglect the Lorentz force as concerns the transversal 
components of the spatial motion. In this respect, we treat the particle as neutral — in 
fact, neutrons would even provide a more natural physical realization of our model than 
electrons. Generally, this is justified as a good approximation if the longitudinal velocity 
is sufficiently small. As an unwanted side effect, with the Lorentz force we also loose 
the time arrow inherent in the interaction with the magnetic field, as concerns orbital 
motion. However, it remains effective through the spin dynamics, see 12.31 below. 

We can therefore restrict the set of relevant dynamical variables to the longitudinal 
coordinate x and its conjugate momentum p, as well as the spin vector s = (s x , s y , s z ). 
In terms of these variables, the Hamiltonian reads 

#(p,x,s;t) = -^- 7 B(x,t)-s, (1) 
zmo 

where mo is the particle's mass and 7 the gyromagnetic ratio. 
We choose the space and time dependence of the field as 

00 

B(x, t) = (0, B 1 {x), B 2 {x)) 6 ( nT ~ *)> ( 2 ) 

n=— 00 

with 

B 1 (x)=A 1 f{x + a/2), B 2 {x)=A 2 f(x-a/2), (3) 
and an envelope function 

e[a/2-|x|], (4) 

as depicted in figure [Jji. This function vanishes outside the interval [—a/2, a/2] yet is 
infinitely often differentiate. 



f(x) = exp 



(a/2? 



or 



2.2. Equations of motion 

In the Hamiltonian (pQ), the spin s plays an exceptional role in that kinematically it 
should be treated as an angular momentum, while there is no corresponding kinetic 
energy term, and a canonically conjugate angle is not defined. For a correct handling of 
these subtleties, we follow the development in [28] and obtain the equations of motion 

m x = 7V(B • s) , s = 7 s x B. (5) 

In all that follows, itlq = 7 = r = 1 is understood. An important feature of (jSJ) is that 
the spin dynamics preserves the form of a Lorentz force and thus breaks time-reversal 
invariance (TRI). 
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The field configuration defined in (jSED implies that within the intervals [—a, 0] and 
[0, a] , the field varies only in magnitude but not in direction. This allows for further 
simplifications. We identify the ^-coordinate with the local field direction, e y = B/fi, 
B = |B|, choose spherical coordinates for the spin, 

s x = |s| sin#siny9, s y = |s| cos 9, s z = |s| sin 9 cos ip, (6) 

and denote s = |s| = 1. Then the ^/-component of the spin, m = s cos 9, is the angular 
momentum canonically conjugate to the azimuth (p, while the polar angle 9 becomes a 
constant of motion, so that the spin dynamics reduces to a mere precession, see figure 

x = mdB/dx <fi = —B (7) 

with m = const. Note that the sense of precession is inverted under time reversal t — > — t 
(keeping B — >• B), revealing the time arrow inherent in ([5]). 




Figure 1. Magnetic field and spin motion in the scattering region, (a) Configuration 
of the magnetic field (|3I4[) . In each of the two sectors 1 and 2, the field is isotropic, 
with an angle of ir/2 between the two sectors. For equal amplitudes and identical 
envelopes, the field is symmetric with respect to rotation by n about the bold line, 
corresponding to the transformation ([9|). Rendering the widths and/or the amplitudes 
of the two field sectors different breaks this symmetry, (b) Schematic representation 
of the scattering process involving the spin vector s (bold arrows). In the incoming 
and outgoing asymptotic regions (phases I, V, resp.) s remains constant. Within the 
field sectors 1 and 2 (phases II and IV, resp.) it precesses around B. Upon passing 
from sector 1 to 2 (phase III), the angles 6, (j> defining the orientation of s undergo a 
passive transformation, cf. ([8|), as depicted in the enlargement (inset). 
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2.3. Symmetry considerations 

For an inhomogeneous but unidirectional magnetic field, according to (JZJ), the projection 
of the spin onto the field direction is a cyclic variable. This kind of dynamics therefore 
may give rise to chaotic spatial motion, but will not affect the spin. In order to induce 
"spin flips" (again abusing quantum terminology), we add a second interaction sector 
where the field is perpendicular to that in the first, as specified in (J3J) (see figured^). 
Upon passing from region 1 to 2 or back, the spherical coordinates (EI) defining the spin 
orientation in the two respective systems of reference aligned with the field direction in 
either region undergo a passive transformation as follows, cf. the inset in figure [lb, 

6*2 = arccos(sin 6*i cos^i), Q\ = arccos(sin 9 2 sin <p 2 ), , > 

y? 2 = arctan(cot Q\ esc <y?i), (pi = arctan(tan# 2 cos<^ 2 )- 

If the widths and amplitudes of the two field sections are identical, however, a spatial 
symmetry remains, namely rotation about the line through the origin and diagonal 
between the two field directions, see figure [Eh- In Cartesian coordinates, it corresponds 
to the operation 

x' = -x, y' = z, z' = y. (9) 

This symmetry would impede directed transport. In order to break it, we choose the 
two amplitudes distinct, Ai ^ A 2 , which gives the difference A 2 — Ax the meaning of a 
symmetry-breaking parameter. 

Besides spatial symmetries, time-reversal invariance deserves some special 
consideration. As discussed in [3ll [32j [33], [34] , it must also be broken to achieve 
directed transport because otherwise the existence of pairs of symmetry-related but 
counterpropagating trajectories leads to an exact cancellation of currents. In [3TJ [32j 
[331 EH [35] , this has been achieved by imposing an asymmetric profile on the driving 
force. In our case, the driving is symmetric under time reversal, while TRI is lifted by the 
presence of a magnetic field. This effect, however, is not obtained through the Lorentz 
force — which indeed we neglect, cf. 12.11 Rather, it is the spin dynamics s = 7 s x B, cf. 
(JSi), that is altered under t — > —t. 



2-4- Discrete-time dynamics 

The impulsive driving allows us to integrate the equations of motion from one kick to 
the next. Placing time sections immediately before each kick, i.e., t n = nr — + , we 
arrive at the following stroboscopic maps for the two respective field regions 

Pn+i = Pn - 2xi, n Bi(xi >n ) cos6 l ii „/[(a/2) 2 - x-J 2 , (10a) 
V?j,„+1 = tpi, n - Bi(x it n) , (105) 

x i,n+l = x i,n + Pn+1 > (10 C ) 

where i = 1,2, refers to (j3J), y^,n refers to the frames related by (JED, and 

%i,n = x n — {— ) % cl/2. One readily checks that these maps are canonical, that is, the 
determinant of their stability matrix equals unity. 
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Figure 2. The time lag t- ln of a scattering trajectory (bold zigzag curve) relative to 
the periodic kicks of the driving (vertical lines) is defined with respect to the moment 
when it enters the scattering region at x = ±a. It can alternatively be scaled as a 
phase shift <j)- m = 27rfi n /r or a spatial offset x m — Pi n t m . By contrast, the definition of 
the delay time refers to the extrapolation of the incoming and outgoing asymptotes 
(dotted lines) forward and backward, resp., till they intersect the line x = 0. See text 
for details. 

In this periodically driven scattering system, the phase shift of an incoming 
trajectory relative to the driving field appears as an additional scattering parameter 
[26l [27] [35] . It can be defined in various manners; for our particular case of an impulsive 
driving and a compact spatial support of the field we choose the time lag of the trajectory 
when it enters the field region at x = ±a with respect to the most recent kick at t — nr, 
cf. figure [2J By definition, its range is < t- m < r. Equivalently, it can be rescaled as a 
spatial shift x in = Pi n t in , < x in < Pi n r, or a phase angle (p in = 2nt in /T, < in < 2tt. 

2.5. Time scales 

The dynamics generated by (JSJ) is characterized by three time scales: (i) the temporal 
separation of the kicks, which defines our unit of time, (ii) the period of spin precession, 
given by the inverse field strength -B -1 , and (Hi) the typical time to pass the scattering 
region, of the order of a/p in . As it is our objective to study the participation of the 
internal freedom in chaotic scattering, we shall work in a regime where all the three 
time scales are comparable. This prevents, in particular, using adiabatic approximations 
based on a slow motion of the external coordinate as compared to spin precession. 

By contrast, in figure [T2b below, we consider the case of fast spin precession 
Ai ^> 1, i — 1,2, in order to anticipate the regime of adiabatic spatial motion typical 
for experiments with spin-| systems. 
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3. Chaotic scattering 

We here consider chaotic scattering as defined by the following properties [lO, SU H2] : (0 
the existence of a chaotic repeller consisting of a discrete set of unstable periodic orbits 
in the scattering region, (ii) rapidly fluctuating deflection functions with self-similar 
structure, at least in a statistical sense, and singularities accumulating towards the 
orbits of the chaotic repeller, and (Hi) an exponential distribution of delay times inside 
the scattering region. In the subsequent paragraphs we present numerical evidence that 
the system devised above complies with all these criteria. 

3.1. Deflection functions 

In plotting the deflection functions (figure [3]) we concentrate on two parameters as 
incoming variables: the phase 0; n (left column, panels a,c,e) as an important control 
parameter in experimental applications and the polar angle 8 m (right column, panels 
b,d,f ) as we are interested in the directed transport of angular momentum and ought to 
make sure that this variable participates in the chaotic scattering. 

While these figures show the typical behavior of a chaotic scatterer for large parts 
of the parameter space, we observe a conspicuously regular pattern in the right column 
of figure [3] (panels b,d,f) for #i n ^7r/2 that calls for a special explanation. It reflects 
a strong asymmetry of the scattering process: Trajectories entering region 1 from the 
left feel a strong repulsive force and bounce back immediately without ever passing 
into region 2, while trajectories coming in from the right undergo the typical irregular 
scattering with unbounded delay time. Such long trajectories tend to randomize the 
outgoing direction and thus lead to approximately balanced transmission and reflection 
probabilities. We shall see in section H] that this asymmetry is largely responsible for 
the transport processes in the system. 

3.2. Stable and unstable periodic orbits 

From fl 1 qj) it is evident that the polar angle 9 plays a decisive role for the dynamics as 
its sign determines the local stability of a trajectory. Specifically, upon passing from one 
field region to the other, the dynamics turns elliptic or hyperbolic according to whether 
cos#in is positive or negative, respectively. This is clearly visible also in figure [6b as an 
absence of peaks for 9 in > tt/2. 

In figure [5] we depict a stable periodic orbit (a) that is inaccessible from outside the 
scattering region and an unstable one (b) that belongs to the chaotic repeller. They are 
marked by symbols + and x, resp., in figure H^. While for the former, the spin remains 
at rest, for the latter it exhibits a significant motion that cannot even be reduced to mere 
precession. In orbital phase space, this orbit forms a double loop, but the corresponding 
self-crossing is lifted upon taking the spin motion into account. 
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Figure 3. Deflection functions. Outgoing azimuth (p out (a,b), outgoing linear 
momentum p out (c,d), and outgoing polar angle 8 out (e,f) vs. initial phase <p- m (a,c,e) 
and vs. initial polar angle 9 m (b,d,f), respectively. The other initial conditions and 
parameters are p m = 1, 9 m = tt/A (a,c,e), p in = 0.5, </>j n = (b,d,f), and <^; n = 0, 
At = 2, A 2 = 1, a = 4. 
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Figure 4. Poincare sections in the (x,p)-plane for positive (panel a) vs. negative (b) 
incoming polar angle 6 ln = — it /A (b). The other initial conditions and parameters 
are (f>- m = 0, tp- m = 0, A\ = A% = 1 and a — 4. The large regular areas in panel a 
correspond to stable islands not accessible from outside the scattering region. Symbols 
+ and x indicate periodic orbits depicted in figure [5k,,b, resp. 



a b 




-10 4 



Figure 5. Examples of a stable (a) and an unstable (b) periodic orbit, located 
respectively in the left stable island in figure 2] (marked by +) and the chaotic region 
(x). The spatial motion is shown as position x vs. continuous time t, the spin 
orientation at each t n in an independent 3-dim. coordinate system moving with the 
spatial trajectory (inset). While the stable orbit leaves the spin at rest parallel to the 
field, the unstable one involves a significant spin motion. 



3.3. Delay-time statistics 

We define the delay time for a given scattering trajectory as usual as the time difference 
between the incoming and outgoing asymptotes, extrapolated forward and backward, 
respectively, till the origin, cf. figure [2j This last specification is important as we are 
dealing with a driven system where in general, energy is not conserved and the outgoing 
momentum is different from the incoming one. 
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Figure 6. Delay-time statistics as a function of initial phase <jy m (a) and incoming 
polar angle 9 ln (b). The other initial conditions and parameters arepi n = 1.0, #i n = 7r/4 
(a), p in = 0.5, <f) in = (b), and ip in = 0, A x = 2, A 2 = 1, a = 4. 



4. Directed transport 

We have shown in the preceding sections that our model fulfills the two necessary 
conditions for directed transport mentioned above, inhomogeneous phase space and 
absence of binary spatio-temporal symmetries, and therefore expect to find directed 
currents at least in parts of its parameter space. 



4-1- Asymmetric scattering and directed currents 

In the context of a scattering system, we define the current as the frequency of particles 
leaving the scattering region to the right minus the frequency of particles leaving to 
the left. Denoting T a p the fraction of particles transmitted from "channel" a to (3 and 
R aa the fraction of particles reflected from a back into a, with a,j3 = "1" (left) or "r" 
(right), the current is given by [35] 

1 = (T^ + Rn-Rn-Trt). (11) 

In the context of quantum electronics, the spin current is given simply as the 
difference between the currents in the two spin states "up" vs. "down", J s = If — 
with the partial currents defined in turn as in (jl~T|) . As an immediate generalization 
to a larger total spin we would write the spin current as a weighted sum of partial 
currents, with the weights given by the projection of the spin onto some appropriate 
direction of reference. This suggests to define the spin current in a classical context, 
where the angular momentum is continuous, as a weighted integral over partial currents 
parameterized by spherical coordinates, 

I s = — d#sin# / dip cos 9 j (9, ip), (12) 

47r Jo J-n 

with the current density 

j(0, V ) = {j lT {0, ip) + j a (0, ip) - j n {9, cp) - Jvl (9, <p)) , (13) 



Chaotic angular-momentum pump 



12 




71/2 % 7t/2 7t 



Figure 7. Qualitative outcome of scattering, i.e., transmission (black) vs. reflection 
(white) as a function of two initial conditions, #i n and pj n , for a configuration with (a) 
and without (b) the spatial symmetry ([9]) . The other initial conditions and parameters 
are A\ = 1 (a), A\ = 2 (b), and Ai = 1, a = 4, </> in = 0, <^i n = 0. 




7t/2 ?i n72 % 



9 in 9 in 

Figure 8. Effective outgoing spin cos# ou t (white = negative, grey = zero, black = 
positive) as a function of incoming polar angle 9- m and linear momentum p ln , for a 
configuration with (a) and without (b) the spatial symmetry ([9]). The other initial 
conditions and parameters are Ai = 1 (a), A\ = 2 (b), and A2 = 1, a = 4, <j5i n = 0, 

<An = 0. 

using obvious shorthands for the partial currents between the respective leads. The 
angles should refer to some suitable laboratory frame, here chosen as the system 
corresponding to field region 1, cf. ([3]). 

We here assume an initially unpolarized ensemble, corresponding to a homogeneous 
angular-momentum distribution over the unit sphere. In some cases it is preferable to 
fix part or all of the initial conditions and to analyze the transport properties as a 
function of these variables. In figures [7] and El we show respectively the outcome of the 
scattering process, i.e., transmission vs. reflection, and the polarization of the outgoing 
particles as functions of the initial polar angle and linear momentum. Also in these 
quantities we observe fractal self-similar structures as in the deflection functions, and 
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a correspondingly sensitive parameter dependence. The presence (panels a) or absence 
(panels b) of the symmetry (Q is clearly reflected in the transport features. The large 
almost void regions that appear in all graphs for 9- m > tt/2 correspond to trajectories 
that bounce back immediately after entry, while manifestly asymmetric structures are 
found only for 6 m < n/2. This confirms the mechanism described in 13.11 that chaotic 
randomization of the outgoing direction occurs only for particles spinning in one sense 
but not for the other. 

The occurrence of directed transport in this system manifestly violates linear- 
response theory. In terms of static quantities, we are dealing with a non-zero mean 
(direct) current in the absence of any external mean gradient. Even if we allow for time- 
dependent input and output, the fact that the system is driven periodically through a 
single parameter would entail an exact cancellation of transport in a linear-response 
treatment. The same conclusion results from a Taylor expansion of the current in terms 
of the driving frequency u, 



If only the ij-term were present (linear response), time reversal of the driving u — ► —uj 
would lead to an exact reversal of the current ("turning the crank backwards pumps 
in the opposite direction"). In our case, the driving is inherently symmetric in time so 
that the presence of a directed current is incompatible with linear response [35]. We 
emphasize that the functional dependence of the current output on the driving amplitude 
is irrelevant for this conclusion. 

The dependence of the polarization on width and amplitude of the field, figure [9h, 
demonstrates how the spin output could be controlled varying these two parameters that 
are easily accessible in experiments. The role of the symmetry is more directly revealed 
in figure [9b where the diagonal A\ = A 2 coincides with a zero of the polarization 
(corresponding to a medium grey shade in the figure). There are, however, other 
conspicuous zones of vanishing polarization apparently unrelated to this symmetry. 

4-2. Angular-momentum transport 

The data presented in figures [8] and [9] strongly indicate the existence of polarized currents 
at least in certain regions of the parameter space. We sketch in figure [TOl how a separation 
of spins could come about in principle in the absence of a net mass current. In order 
to corroborate its robustness and experimental feasibility, we show in figures [11], [12] 
global averages both of the particle and the spin current, cf. flll |12p . We find zeros 
of the spin current I s at appreciable values of the particle current J, corresponding to 
unpolarized charge transport, as well as zeros of I at high values of I s . These latter 
cases amount to a separation of different orientations of angular momentum without net 
particle transport, that is, to the classical analogue of spin filtering or rectification. 

Moreover, in figure [H] we observe a strong dependence of this phenomenon on the 
breaking of parity, in terms of the ratio A\/A 2 - As discussed in 12.31 above, it must 
vanish for A\ — A%. At the same time, we expect it to diminish in the opposite extreme 




(14) 
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Figure 9. Effective outgoing spin cos6* ut (white = negative, grey = zero, black = 
positive) as a function of width a and amplitude Ax of the left field sector (a) and of 
the amplitudes Ai and A2 (b), averaged over 0- ln and p m - The other initial conditions 
and parameters are A 2 — 1 (a), a = 4 (b), 0i n = and tp- ln = 0. 




Figure 10. A finite spin current in absence of net mass (charge) transport can 
arise if particle currents cancel due to parity invariance between transmission and 
reflection from either side, while the corresponding coefficients for spin scattering lack 
this symmetry. In this schematic figure, all particles with spin down (full lines) leave 
to the left, irrespective of the incoming direction, while those with spin up (dotted) all 
leave towards the right. At the same time, if all currents are assumed to have the same 
magnitude, mass transport vanishes identically, resulting in pure spin separation. 



of Ai <C A 2 or v.v., as in this limit the effect of the two regions with different field 
orientation breaking the symmetry is lost. Therefore there should exist an optimum 
for angular momentum separation at intermediate values of Ai/A 2 . 

5. Conclusion 

With this work we intend to demonstrate the feasibility of pumping angular momentum 
in a chaotic scattering system, in the regime of fast and strong driving where adiabatic 
or perturbative methods would fail. The transport phenomena we observe, both of mass 
and angular momentum, owe themselves to an interplay of strongly nonlinear dynamics 
and the breaking of spatio-temporal symmetries. They go far beyond the frame of 
linear response and do not require a two-parameter driving as in adiabatic pumping. As 
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Figure 11. Particle current / (panel a) and spin current I s (b), averaged over 9 m and 
Pi n vs. the width of the well a for various values of the ratio Ai/A 2 . Parameters are 
A 2 = 1 and ipi n = 0. 





Figure 12. Comparison of particle current / (broken curve) to spin current I s (dotted), 
averaged over 9 ln and p lni as functions of the width a of the scattering region (panel a) 
and the mean magnetic field strength (panel b). Parameters are ip- la = and A\ = 8, 
A 2 = 1 (panel a) and a = 4, Ai/A 2 = 1.5 (b). 



a welcome side effect, the sensitive parameter dependence of the scattering process opens 
the possibility not only to generate angular-momentum currents but also to control them 
on short time scales and thus to involve them in computation processes. 

At the same time, our results constitute an important example of how an internal 
degree of freedom not only participates in chaotic scattering, but even gives rise to new 
phenomena. While in this work we focus on transport, the aspect of chaotic scattering 
with internal freedoms deserves a closer scrutiny in a separate effort. 

Another relevant side issue is the fact that spin separation clearly corresponds 
to a reduction of entropy. In this paper we remain within a perfectly deterministic 
Hamiltonian setup, a more realistic model including thermodynamic aspects however 
would urgently have to address this question. 
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Similarly, connecting the pump on both sides with reservoirs as is usually considered 
in the context of electronic transport phenomena opens additional possibilities to 
achieve directed currents. The distributions characterizing the reservoirs (e.g., Fermi 
distributions parameterized by the chemical potential) require dissipative processes to 
remain invariant and thus implicitly constitute a time arrow. As a consequence, in such 
a configuration directed transport can arise even if the scattering system itself does not 
break time- reversal invariance. 

The classical study presented here can be considered as a preparatory survey aiming 
towards chaotic spin pumping. On grounds of semiclassical arguments, we would expect 
the mechanisms we elucidate to carry over at least to a "mildly" quantum regime. 
The electron spin, however, corresponds to the deep quantum limit where among other 
consequences, the assumption of a slow spin precession, comparable to the other time 
scales in the system, no longer applies. In order to address this question, we depict in 
figure [12b mass and spin currents as a function of the magnetic field strength, equivalent 
to the precession frequency, at constant asymmetry. We see that directed transport does 
not diminish significantly and we even find instances of pure spin current in a regime 
where both time scales differ by more than an order of magnitude, giving us some 
confidence that chaotic pumping of angular momentum could extend down to s = H/2. 
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Abstract. Chaotic scattering with an internal degree of freedom and the possibility 
to generate directed transport of angular momentum is studied in a specific model, 
a magnetic dipole moving in a periodically modulated magnetic field confined to 
a compact region in space. We show that this system is an irregular scatterer in 
large parts of its parameter space. If in addition all spatio-temporal symmetries are 
broken, directed transport of mass as well as angular momentum occurs. The sensitive 
parameter dependence of the corresponding currents includes frequent sign reversals. 
Zeros of either quantity entail the exclusive occurrence of the other and thus give rise 
in particular to angular-momentum separation without mass transport as a classical 
analogue of spin-polarized currents. 
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1. Introduction 

During two decades of research on chaotic scattering, internal degrees of freedom have 
not enjoyed much interest. Yet there are numerous reasons to give them a closer look: An 
internal freedom may provide the additional dimension required to render a scattering 
system with a single external degree of freedom chaotic [?, ?, ?, ?]. It may also act as a 
reservoir that absorbs energy from the external motion and stores it temporarily thus 
giving the scattering process a transient inelastic character. Concerning applications, it 
is obvious that in chemical reactions — one of the paradigms of chaotic scattering [?, ?]- 
the presence of internal freedoms is indispensable. In mesoscopic physics, in turn, the 
basic entities are electrons or fermionic quasi-particles: Their spin gives rise to a rich 
phenomenology beyond point-particle physics and presently receives much attention in 
the context of spintronics [?, ?, ?]. 

In electronic and other mesoscopic systems, transport is an issue of vital importance. 
It has been studied in depth in the framework of adiabatic pumps [?, ?, ?, ?], electron 
counting [?, ?] and of the Kubo and Buttiker formulae [?, ?, ?]. In spintronics, the 
practical necessity to provide polarized currents has spurred research on directed spin 
transport [?,?,?]. In a different context, the study of ratchets [?], inspired by the 
biological phenomenon of motor molecules [?, ?, ?], has revealed general physical 
conditions and numerous mechanisms for directed transport to occur in systems far 
from thermal equilibrium. Also here, internal freedoms are almost always involved 
and participate actively in transport processes [?], e.g., as functional (conflgurational) 
degrees of freedom in motor molecules [?, ?] and as internal heat baths. 

In this paper we present a case study of chaotic scattering with an internal degree 
of freedom, intrinsic angular momentum, and consider how it contributes to irregular 
scattering and how in turn directed transport of this quantity arises from scattering. 
The subject bears on various of the research lines mentioned above: We shall invoke 
much of what is known about chaotic scattering, in particular in periodically driven 
systems [?, ?]. The coupling of a magnetic dipole to an inhomogeneous magnetic field 
[?] is crucial for our modelling. Also concerning chaotic dynamics in such a configuration 
we refer to earlier work [?]; in fact, it is the essential ingredient of one of the first models 
proposed for irregular scattering [?]. 

In directed transport, we take all the insights into account that concern the 
importance of binary spatio-temporal symmetries and their breaking [?, ?, ?, ?]. We 
can here fall back on direct precursor work on chaotic pumps [?]. The definition of 
spin current is adopted from quantum contexts to define angular-momentum transport 
in a classical setting. It should be emphasized that we consider a driving that 
is both strong and fast [?, ?] and thus not amenable to any perturbative and/or 
adiabatic approximation. In this sense, we are dealing with strongly nonlinear transport 
phenomena, far off the regime of linear response and in particular of peristaltic pumping 
[?, ?, ?, ?]. As a consequence, directed transport is achieved already driving the system 
via a single parameter. The crucial role of chaos also distinguishes our work from recent 
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studies of spin ratchets in ballistic electron systems [?, ?, ?]. Preliminary results have 
been published on a national platform [?]. 

As our main achievement, we present conclusive evidence that in classical chaotic 
scattering systems, polarized currents can be generated, i.e., directed transport of 
angular momentum occurs independently of and even in the absence of mass (charge) 
transport. We explain this effect in terms of the interplay of asymmetry and the 
randomizing action of long unstable trajectories, thus corroborating the decisive role 
of chaotic scattering. It is also manifest in the sensitive parameter dependence of mass 
as well as angular-momentum currents, which allows for a fine tuning of both quantities. 
In particular, sign reversals are frequent and zeros of one current give rise to the exclusive 
dominance of the other. In this way, we achieve a "rectification of angular momentum" 
in the absence of mass transport. 

In ??, we construct our model and justify its peculiarities. We discuss details 
of the dipole-field coupling that have not been considered previously in the context 
of chaotic scattering. Section ?? is dedicated to irregular scattering as the basic 
dynamical category of the system. Various diagnostics are presented such as deflection 
functions, unstable periodic orbits, and time-delay statistics. Our central subject, 
directed transport, is approached in ??. We define the relevant transport quantities 
in terms of reflection and transmission coefficients. Numerical evidence is provided 
for directed mass as well as for angular- momentum transport, and we point out the 
manifolds in parameter space where pure currents of either kind occur. Section ?? 
contains a summary and an outlook towards the quantization of our model with the 
perspective to construct a chaotic spin pump — evidently a major motivation for the 
present purely classical study. 

2. A model for chaotic scattering with intrinsic angular momentum 

2.1. Basic setup — coupling internal to external freedoms 

In constructing our model, we attempt to satisfy the following requirements, that (i) 
the system show chaotic scattering at least in parts of its phase space, as a dynamical 
mechanism for directed transport, (ii) the potential couple the internal to the external 
degree of freedom and (Hi) vanish identically outside a compact scattering region so that 
numerical simulations can be restricted to a spatial box of finite extension, and (iv) be 
periodically time dependent in the form of delta kicks in order to facilitate reducing the 
dynamics to discrete time, that is, a map. 

Moreover, previous work on nonlinear transport mechanisms in ratchets and pumps 
has revealed a fundamental condition for directed currents to occur [?, ?, ?, ?, ?]: an 
inhomogeneous phase space reflecting the absence of any spatio-temporal symmetry that 
would give rise to pairs of otherwise identical trajectories carrying the same current 
in opposite directions. It can be accomplished in particular if the dynamic is mixed 
(regular coexisting with chaotic motion) and the associated invariant manifolds exhibit 
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some asymmetry. We take these additional criterion into account in setting up the 
model. 

Specifically, we choose the internal freedom as an intrinsic angular momentum (in 
the following we shall also refer to it as spin, wherever its classical nature is obvious). 
In connection with a charge it gives rise to a magnetic dipole moment that couples 
to an external magnetic field and thus, if the field is not homogeneous, to the spatial 
motion of the particle. However, we neglect the Lorentz force as concerns the transversal 
components of the spatial motion. In this respect, we treat the particle as neutral — in 
fact, neutrons would even provide a more natural physical realization of our model than 
electrons. Generally, this is justified as a good approximation if the longitudinal velocity 
is sufficiently small. As an unwanted side effect, with the Lorentz force we also loose 
the time arrow inherent in the interaction with the magnetic field, as concerns orbital 
motion. However, it remains effective through the spin dynamics, see ?? below. 

We can therefore restrict the set of relevant dynamical variables to the longitudinal 
coordinate x and its conjugate momentum p, as well as the spin vector s = (s x , s y , s z ). 
In terms of these variables, the Hamiltonian reads 

H(p,x,s;t) = 7 f-->yB(x,t)-s, (1) 

where mo is the particle's mass and 7 the gyromagnetic ratio. 
We choose the space and time dependence of the field as 

00 

B(x,t) = (0,B 1 (x),B 2 (x)) S(nr-t), (2) 

n=— oo 

with 

B x (x) = AJ(x + a/2), B 2 (x) = A 2 f(x - a/2), (3) 
and an envelope function 

9 [a/2 - |x|] , (4) 

as depicted in figure ??a. This function vanishes outside the interval [—a/2, a/2] yet is 
infinitely often differentiable. 



f(x) = exp 



(«/2) 5 
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2.2. Equations of motion 

In the Hamiltonian (??), the spin s plays an exceptional role in that kinematically it 
should be treated as an angular momentum, while there is no corresponding kinetic 
energy term, and a canonically conjugate angle is not defined. For a correct handling 
of these subtleties, we follow the development in [?] and obtain the equations of motion 

m x = 7V(B • s) , s = 7 s x B. (5) 

In all that follows, m = 7 = r = 1 is understood. An important feature of (??) is that 
the spin dynamics preserves the form of a Lorentz force and thus breaks time-reversal 
invariance (TRI). 



